function obj = A1a_mc_f_opt_dom(mc_f, EstShare_true, price_k_market, expenditure_k_market, alpha_price, alpha_1, Y_market, simShare_true, nn, N, pos_dom)

mc = mc_f(1:(end-1),:);
f = mc_f(end,:);

N_j = EstShare_true.*(price_k_market - mc) ./ price_k_market .* expenditure_k_market/f;
alpha_i = alpha_price + alpha_1*Y_market;
sum_alpha_s = alpha_i.*simShare_true;
sum_alpha_s2 = alpha_i.*(simShare_true.^2);
obj1 = EstShare_true./N_j .* mc + (price_k_market - mc) .* (1/nn) .* [sum(sum_alpha_s,2)./N_j - sum(sum_alpha_s2,2)./(N_j.^2)];   
obj2 = N - sum(EstShare_true(pos_dom, 1) .* (price_k_market(pos_dom, 1) - mc(pos_dom, 1)) ./ price_k_market(pos_dom, 1) .* expenditure_k_market/f);
obj = [obj1; obj2];

end

